Role of the Landau-Migdal Parameters with the Pseudovector and the Tensor 

Coupling in Relativistic Nuclear Models 
— The Quenching of the Gamow- Teller Strength — 
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Role of the Landau-Migdal parameters with the pseudovector (g a ) and the tensor coupling (gt) is 
examined for the giant Gamow- Teller (GT) states in the relativistic random phase approximation 
(RPA). The excitation energy is dominated by both g a and g t in a similar way, while the GT strength 
is independent of g a and gt in the RPA of the nucleon space, and is quenched, compared with that 
in non-relativistic one. The coupling of the particle-hole states with nucleon-antinucleon states is 
expected to quench the GT strength further through g a . 

PACS numbers: 21.60.Jz, 21.65.+I, 24.30.Gd 



I. INTRODUCTION 



The long- and medium-range parts of the nuclear spin- 
isospin interaction are well described by the one-meson 
exchange approximation. In the non-relativistic limit, 
they are written in the momentum space asQ 
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Jt L \ (o-i • q)(o-2 • q) 
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where the first term stands for the one-pion exchange 
potential, while the second term the one-rho exchange 
potential. They are obtained from the pseudovector 
and tensor meson-nucleon couplings, respectively. In the 
short-range part of the interaction, on the other hand, 
many-body correlations should be taken into account. 
Their effects are approximately expressed by the Landau- 
Migdal (LM) parameter g' 0, 
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Experimentally the value of g' is estimated to be about 
0.60 ■ This fact is understood in such a way that about 
a half of the value cancels the zero-range part of the one- 
pion exchange potential, while the rest the short-range 
part of the rho-meson exchange potential in Eq. 0, |4j • 
The contribution from the exchange terms of the one-pion 
and one-rho potential is also considered to be included 
in g' . Indeed, if we calculate the LM parameter so as to 
cancel the zero-range part of the one-pion potential and 
to include the contribution from the exchange term of 
the finite-range part, we obtain from Eq.Q 
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This value is about 0.318 for the pion mass m„ = 140 
MeV and the Fermi momentum fep = 1.36fm~ . The 
corresponding expression from the rho-meson exchange 
potential is given by 
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If we simply calculate this value by using the strong cou- 
pling / 2 /4tt = 4.86 together with m p = 770 MeV, we 
obtain g' p = 1.07. Refs.Q and 0] have shown, however, 
that in the p-meson exchange, the short-range part of the 
finite range interaction should be also taken out, since the 
p-mass is much larger than the 7r-mass. This contribu- 
tion makes the above value of g' p smaller as 0.281, since 
the finite range part has an opposite sign to that of the 
zero-range part. Thus the experimental value about 0.6 
is well understood by the sum of g'^+g' in non-relativistic 
models. 

Now, for the last 30 years it has been shown that nu- 
clear structure is well explained phenomcnologically by 
relativistic modelsQ. Not only the ground state, but also 
the excited states are well reproduced ,6j. So far, however, 
spin-dependent excitations of nuclei have not been stud- 
ied in details 0. One of the reasons may be because there 
is no estimation of the LM parameter yet in relativistic 
models. Moreover, there is no unique way to extend the 
above non-relativistic g' to the relativistic one. 

The purpose of the present paper is to introduce phe- 
nomenologically the LM parameter in relativistic models, 
and to explore model-dependence of the excitation energy 
and strength of the giant Gamow- Teller (GT) resonance 
state. We will construct relativistic LM parameters so 
as to reproduce g' of Eq.J2J) in the non-relativistic limit. 
Such candidates will be given by contact terms with the 
pseudovector (g a ) and the tensor (g t ) coupling in the nu- 
clear Lagrangian. 

In the previous paper 0, 0, El, we investigated al- 
ready the LM parameter of the pseudovector type for 
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the description of the giant GT resonance states. It was 
shown that the excitation energy in non-relativistic mod- 
els is reproduced and that the GT excitation strength 
in the nucleon space is quenched, compared with the 
sum rule value by Ikeda-Fujii-Fujita (IFF)^]. In par- 
ticular, the quenching, which is consistent with recent 
experiment 12], was an important conclusion, since in 
conventional non-relativistic models the value of the GT 
strength is fixed by the IFF sum rule. In this paper we 
will pay our attention mainly to whether or not those 
conclusions in the relativistic model are changed by the 
LM parameter with the tensor coupling. 

In the following section, we will calculate the correla- 
tion function in the random phase approximation (RPA), 
where the correlations are induced through the LM pa- 
rameters g a and g t . In Sec lIIII and IIVI the excitation 
energy and strength of the GT state will be estimated, 
respectively. Since all calculations are performed for nu- 
clear matter, discussions will be given according to an- 
alytic formulae of the excitation energy and strength. 
They make clear structure of the giant GT states in rel- 
ativistic models. The final section will be devoted to a 
brief conclusion of the present paper. 



II. RELATIVISTIC RPA CORRELATION 
FUNCTION 

We assume that nucleons and antinucleons in nuclear 
matter are bounded in the mean field of the Lorentz 
scalar and vector potentials. The Lorentz scalar poten- 
tial will be included into the nucleon effective mass M* 
below. On the other hand, we will not describe explicitly 
the Lorentz vector potential throughout the present pa- 
per, since it does not play any role in description of the 
GT states. 

The particle-hole and nucleon-antinucleon correlations 
are assumed to be caused by the Lagrangian, 



(5) 



where the first term describes the pseudovector coupling 
and the second one the tensor coupling, 

We have also defined the isospin operators for conve- 
nience as 



t x ± ir, 
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TO = T z 



The coupling constants g a and gt will be related later to 
the LM parameter g' in non-relativistic models. 

The RPA correlation function for the external field, r A 
and Fb, is given by the mean field correlation function 
II as 
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n RPA {r A ,r B ) = n(r A ,r B ) 

+ Xa n(r A , T? ) n RPA {r m ,r B ), (6) 



where we have used the notations, 
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(a = a) 
,, (a = t). 



When the external fields are written as r A = ^aJ- , an d 
r B = 7b t+, 7a and 7b being some 4x4 matrixes, Eq.JBJ) 
becomes to be 

#RPA(7aT-,7bT+) 

= 77(7 a T„,7 b T + ) 

+ Xai7(7 a T_,r^)77 R p A (/V,7 b T + ). (7) 

In the same way, the last term of the above equation is 
described as, 

^?RPA(^Q-,7bT+) 

= i7(r Q _, 7b r+) 
+xpU{r a -,if )i7 R p A (r /3 _, 7b T + ). (8) 

This can be rewritten as, 

iTRPA0Ta-,7bT + ) - (t/- 1 ) ct/3 i7(r /3 _, 7b T + ), (9) 
with use of the dimesic function, 

u afi = 5 a p-xpn(r a _,rl). (10) 

From Eq.©, Eq.Q is expressed in the form, 

^RPA(7aT_,7bT+) 
= 77(7 a T_,7bT + ) 

+ X ai7(7aT-,r«) (C/- 1 ) Q/3 7T(r /3 _,7bT + ). (11) 

Thus, TJrpa can be described in terms of the mean field 
correlation function 77. 

When we calculate the mean field correlation func- 
tion, we neglect the divergent terms and keep all the nu- 
clear density dependent terms including the Pauli block- 
ing ones, as most of the previous authors did[L3l fi~4| . 
We call this approximation no free term approximation 
(NFA)[l(J. As discussed in details in Refold], the Pauli 
blocking terms are necessary for keeping at least the IFF 
sum rule and the current conservation|lOt Il5| . The no- 
sea approximation (NS A) 0] is equivalent to NFA in the 
description of the low lying states where nucleon-degrees 
of freedom play a main role^(j- In the following calcula- 
tions, we will use NFA, but will come back to the problem 
of this approximation at the end of Sec lIVI 

In NFA, the mean field correlation function is given by 

n(-f a T-,-fbT+ ) 

d 3 p ( t ab (p,q) (n) 



(p + qf - M* 2 + is p 



tba(p,-q) 



(p - q) 2 - M* 2 + ie 
3 S(q + E p - £ p+q ) 



0(p) 
p 



+ i«ld*p ~^ ~™' t,6(p,g)gWffi q , (12) 
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where q stands for the four-momentum transfer from the 
external field to the nucleus. The notations are defined as 
Po = E p = \/v 2 + M* 2 , M* being the nucleon effective 



mass, and 6 P = 6(ki — |p|), k n and k p being the Fermi 
momentum of the neutrons and protons, respectively. 
The step functions express the density-dependence. We 
have also used the abbreviation in the above equation, 

tabip, q) = - Tt* (laty + 4+M*) lb y+ M *)) . (13) 

In the present paper, we study the GT transition at q = 
0. In this case, Eq. l|13[l can be written, as shown later 
explicitly, in the following form, 



tabip, q) = fip) +gip)qo- 



(14) 



Then, Eq. H12|) . for N > Z nuclei, can be expressed in 
terms of fip) and gip) as 



ni la T- 

d 3 p 

2 



• lbT+ ) 
f 



+ (2po.9 - /) 



is 



9 (n) 
7 P 



9 (p) 
7 P 



2po + qo - is 2p - q - ie 



For giant GT states with \qo\ -C M*, the above equation 
may be written as 



7T(7 a T_,7 h T + ) = 



d 3 p 

2po L qo + i - 

2poff - / 
2po 



ie V p 



0(p) 



■0$ 
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As shown in the previous paper |9|, the GT states can 
be described by taking the component, j a — 7& = 7572, 
as the external field. This component has non-zero values 
of t a b in the following correlation functions, 

Tin = i7(7572T-,7572T+) 

= -^(7572T-,757 2r +) (16) 

1JX 2 = i7(757 2 T_,CT3lT + ) 

= ^7(7572 r_,(T 31 r + ) 

= - ^(o-3iT_,757 2 r + ). (17) 

Therefore, in writing 

n 22 = n(a 31 T^,a 31 T+) = n(a 31 T^ 1 a 31 r + ), (18) 

the dimesic function is provided as 

u = s aP - X0 n(r a _,rP) 

1 + Xallu ~ XtIIi2 
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1 - XtTl; 



22 



(19) 



Calculations of t a b for 11^ are straightforward. For 
U11, we obtain 

t ab = - 4 (2A/* 2 + 2p 2 + E p q Q ) . 



This, together with Ea. (|15|) . gives 



n n = 



qo + is 



2 d 3 p 



d 3 p 



M* 2 + pi ( , s , s 
tJL [ f)( n ) _ fl(p) 



p2 V P P 

P 



2 2 



(20) 



In the case of 77i2, i a t> is calculated to be 
t o6 = -4M* (2^ + go), 

which yields 



77 



12 



4M* 



3„ °P °P 



d 3 p 



(21) 



In the same way, U22 is described as 



4 f o M* 2 + 2p 2 

qo + ie J 



(,(n) 



by using 



2/^B!_M(,(n) + ,(p)), (22) 



U = - 4 (2M* 2 + 4p 2 + £ P g ) 



In order to express the above equations in a simpler 
form, let us employ the following notations, 



pk F fyj*2 _|_ 2 

?.(**) = 4^ d 3 p^-^L 
M* 



Q m (k F ) = 4 / d d p 



-Bp 



Qt(fcp) = 4 / Vp 
Jo 



E 2 P 



(23) 
(24) 
(25) 



" = ^/ d ^||( p n) +< ) )' (26) 



and, moreover, 



q a = Qaikn) ~ Qa(k p ) , (27) 

(M > (28) 
« = Qt(k n ) - Qtik P ). (29) 



Then, 77^ is described as 
77n = 

n 12 -- 
n 2 2 -- 



q„ 



qo + is 
qm 

qo + i£ 
qt 

qo + is 



2k, 



K. 



(30) 
(31) 
(32) 
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The functions in Ea. (|27() can be expanded in terms of 
(fc n — fc p ), for example, as 

^dg(M _ dQ(M _ 2 / _ 2 2 \ 

9q dfc F (fen * pj ' dfc F - lb7rfc F^l 3 %J, 

where wf stands for the Fermi velocity, fc F / Ep, with Ep = 
yjkl +M* 2 . When we use the relationship, as usual, 

2, N — Z 



k n fcp ps - fcp — , 

the function q a is written in terms of (N — Z) as 

(2 \ N — Z 

1 - -4 J , a = 327rfc|^-. (33) 

In the same way, the functions, q m and qt, are expressed 
approximately as 



'/,„ - n v / I - rp . (/, '-a a ( 1 - ^vp 



(34) 



We note again that the above equations are obtained by 
their expansion in terms of (N — Z), but not in terms 
of vp. On the other hand, if we expand k in Ea. i|26|) in 
terms of vf, we obtain 



8tt 2 o / 3 2 



(35) 



This comes from Pauli blocking terms due to nucleon- 
antinucleon excitations [fl flfTl . Since n is of order tip, the 
Pauli blocking effect is negligible in the case of the GT 
excitations j^, llOf . 



III. EXCITATION ENERGY OF THE GT 
STATES 

The excitation energy of the GT state is determined 
by the determinant of the dimesic function, 

det U = 0, (36) 

which gives go = LU ± ■ 

i (oj a + ut ± V K - ^t) 2 - ^XaXtql, ) (37) 



UJ± 



with 



= Xaqa , Ut = ~ XtQt 

Xa - Xt 

Xa = 7T— > Xt ^ 



1 - 2kXo 



1 + KXt 



(38) 
(39) 



It is obvious that ui a and oj t are the solutions of Eci. (|36[) . 
when there is no mixing between the pseudovector and 
tensor couplings. According to Eas. (|3*3l) and (|3"P . they 
are given approximately as 

u a ~ axa (l - > w t ~ ~ a Xt U - ~4j ■ (40) 



In non-relativistic models, it has been shown that the 
LM parameter a' p rovides us with the excitation energy 
of the GT statejl3, 



4fcp N - Z 
3^ A ' 



(41) 



Comparing the above equation to Eq. l|40|l . it is seen that 
both pseudovector and tensor coupling can reproduce in- 
dividually the non-relativistic result by the relationship, 



9a = g 



9t 9 ) • (42! 



but they produce a different rclativistic correction of or- 
der v F from each other. The above relationship is also 
verified from the comparison of Eq. (J2J) with the one from 
the non-relativistic reduction of the space part of the La- 
grangian Eq.JSJl, 



£ ~ - 9a 2 — ift<rni> ■ ip^erTit/). 



(43) 



If we take both the pseudovector and tensor couplings, 
uj + in Eq. 13 7|) up to second order of v-p is written as 



w+ « a \ Xa ~ Xt 



2Xa - Xt 



while oj_ is of fourth order, 



10- 



%XaXt 4 
:V F . 



(44) 



(45) 



9(Xa - Xt) 

As mentioned in SecQJ the non-relativistic LM parame- 
ter works in a way that a half of g' is for the 7r-meson 
exchange potential, and the rest is for the p-meson ex- 
change one. In this sense, it may be reasonable to use the 
pseudovector and the tensor coupling with equal weight 
in relativistic models. In assuming that 



g a 



-at = - 2 g 



the above uj + becomes 



WGT ( 1 - ~V F ) . 



(46) 



(47) 



IV. THE EXCITATION STRENGTH OF THE 
GT STATES 

The GT strength is given by the imaginary part of the 
RPA correlation function. The RPA correlation function 
is written, using Eq.Q, as 



77, 



RPA 



(7572 T-, 7572 t+) = (U 1 ) lQ iT(r , «_,757 2 r + ). 



The explicit expression of the dimesic function in Eq. I|19|) 
yields 

#rpa(7572T--j7572T+) = -|- fl - d ^ 22 ^j ■ ( 48 ) 
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The last term in the parentheses can be expressed in 
terms of the eigenvalues in the preceding section, 

1 - x±n 2 2 



detC/ 
1 

1 - 2«Xa 



lu + — lu- qo — lu+ + is 
lo + — lu- qq — u>- + is 



(49) 



The response function for the external field, 75 72"? + , 
therefore, is described as 



Im 77 RPA 



Ra(qo) 

3 A 
16^ fcf 

3 A 1 f u) a - 

16n fc F Xa (1 - 2KXa) - 
UJ + -LL> a 



-uJ+5(q Q - W+) 



Lo_5(q -Lu-) . (50) 



The above equation provides us with the GT strengths 
of the two states with the excitation energy, o>+ and to _ , 
respectively, 



3 A 



167T fc F Xa (1 - 2KXa) ^+ - W- 

3 4 1 cij + — cj q 

167T fc| Xa (1 - 2KXa) ^+ - tO- 



LO+ , (51) 



(52) 



If we take into account order up to vp, they become 
S - 3 — ^ a 

+ ~ 167T/c F Xa(l-2KXa) 
1 



(1 - ^XaY 
S- M 0. 



wp ) 2 (JV - Z) , (53) 
(54) 



Thus, the strength of the GT state does not depend on 
the tensor coupling, up to up, although the excitation 
energy does, as shown in the preceding section. The sum 
of the two strengths in Eas. i|51[l and (|52|l is given, without 
expansion in terms of v-p , as 



S- = 



dq R a (qa) 
3 A 



167T kp Xa (1 - 2ftXa) 

1 ( . 2 



(1 - ^XaY 



l--u F )2(N-Z), (55) 



which is independent of the tensor coupling. If we dis- 
cuss the GT strength in the nucleon space only, then the 
sum is independent of g a also, because of k = 0, and is 
quenched by the factor (1 — 2up/3), compared with IFF 



sum rule value 2(N — Z) in the present definition^ ITo| . 
The last equation is what we obtained in our previous 
papergHHjuJ. 

It may be worthwhile noting the response function for 
the external field, 0311"+. The RPA correlation function 
in this case can be written as 

^RPA (°31 T— j "31 T+ ) 
K 



1 + «Xt 

1 f U>t — U3— LJ-i 



Xt (1 + KXt) \u+ - uj_ q - lu + + is 

LO + — LJ t LU- 

LO+ — LU- qo — LO- + is 

This gives the response function, 



(56) 



Rt{qo) 

3 A 
16tt 2 fcf 

3 A 



Im TTrpa 



1 



UJ t - LU- 

16tt fcf xt (1 + K Xt) \w+ - 

LO + - UJ t 



■ u+ S(qo - uj + ) 



UJ+ — US- 



UI- S(q - U-) . (57) 



The sum of the strengths expanded in terms of (N — Z) , 
therefore, is obtained as 



dqo Rt(qo) 



3 A 



u>t 



16tt fc| xt (1 
1 I 



(i+Kxtr 



K Xt) 

1 



which is independent of the pseudovector coupling. As 
expected, the relativistic correction in Eq.|J5HJ| is different 
from that in Eq.JSSJl. 

Finally effects of the Dirac sea in the relativistic model 
should be mentioned in more detail. In the above calcula- 
tions, effects of the Dirac sea are included in n in Ea. (|26[) . 
which comes from the Pauli blocking terms in NFA. Since 
its value is positive, the GT strength in Eas. (|53l) and 
is a little increased owing to the coupling of the 
particle- hole states with the nucleon- antinucleon ones. 
This increasing is, however, due to a poor approximation 
of NFA where the divergent term is simply neglected. 
The no-sea approximation, which has been extensively 
used so for, is essentially the same as NFA, and provides 
us with Eqs.J53J) and as shown in Ref.[T(|. If we 

keep the divergent term in the RPA correlation function, 
k in Eq. (|26[) is replaced by k' [l(| , 



d 3 p 



El 



(59) 



Thus, the divergent term has an opposite sign to that 
of the density-dependent part. This fact implies that if 
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we treat the divergent terms properly, the GT strength 
may be quenched more than in Eqs.JS3J) and l|55l) . In 
this sense, the factor (1 — 2u F /3) yields the minimum 
quenching of the present relativistic model. We note that 
in nuclear matter the quenching due to the antinuclcon- 
degrees of freedom is caused only through the pseudovec- 
tor coupling. In finite nuclei, where the momentum is not 
a good quantum number, there may be a small contribu- 
tion from the tensor coupling. It is important future work 
to investigate how much the quenching is increased by 
the coupling of the particle-hole states with the nucleon- 
antinuclcon states. In that case, it may be also required 
to take into account the energy-dependence of the LM 
parameter which has not been well studied yet. 

V. CONCLUSIONS 

In the previous papers[3,ElEl, we investigated the ex- 
citation energy and strength of the giant Gamow- Teller 
(GT) states in the relativistic model by introducing the 
Landau-Migdal (LM) parameter g a with the pseudovec- 
tor coupling. The pseudovector coupling is chosen so 
as to reproduce the non-relativistic LM parameter g' in 
the non-relativistic limit. The main conclusions were 
that the relativistic correction to the excitation energy 
is 14% and that the GT strength in the nucleon sector is 
quenched by 12% in nuclear matter, compared with the 
Ikeda-Fujii-Fujita (IFF) sum rule value. The quenching 
factor is given by (1 — 2vp/3), vp being the Fermi mo- 
mentum. In finite nuclei, the quenched amount has been 
estimated to be about 6% in Ref . [9l ITsI] . The reduction of 
the quenching is due to the larger effective mass near the 
nuclear surface than in nuclear matter. Among the above 
results, the quenching of the strength is especially im- 
portant. This prediction is a rare example to distinguish 
the relativistic model from conventional non-relativistic 
models, and is related to an important observation of the 
recent experiment [L2| . 

Recently, Tokyo group has observed 10% quenching of 
the IFF sum rule value in 90 Nb through charge- exchange 
reaction|l2j. Since nuclear models assuming the nu- 
cleus to be composed of nucleons satisfy the model- 
independent sum rule, the observed quenching implies 
that additional degrees of freedom play a role to reduce 
the GT strength. So far, one candidate has been pro- 
posed, which is the A-degrees of freedom 0, 0- Particle- 
hole states are assumed to couple with A-hole states 
through the Landau-Migdal (LM) parameter g' NA . Since 
the coupling force is expected to be repulsive, a part of 
the strength in the nucleon space is taken by the highly 
excited states. Of course, it depends on the value of g' NA 
how much strength is taken out from the low excitation 
energy region. 

The relativistic model predicted another source of the 



quenching. If the 6% quenching is due to the relativistic 
effect, as estimated in the previous papers 0, then 
it is about 4% what is left for the coupling with the A- 
hole states. In this case, the value of g' NA becomes less 
than 0.2, which is much smaller than the one estimated 
before 0, [J 113 • The small value changes our previous 
understanding of the spin-dependent structure of nuclei. 
For example, the critical density of the pion condensation 
becomes less than two times of the normal density |20j. 
Thus, it was required to investigate model-dependence 
of the previous result m the relativistic model. 

In the present paper, we have explored whether or not 
the LM parameter with the tensor coupling changes the 
previous conclusions. The tensor coupling is another can- 
didate which can reproduce the non-relativistic result of 
the GT states. We have described the GT states in the 
two cases. One is to use the tensor coupling instead of the 
pseudovector coupling. The other is to take into account 
both pseudovector and tensor couplings. 

If we use the only tensor coupling, the relativistic effect 
on the excitation energy is reduced by a half, as shown in 
Eq. ■ Eq. iJSBJ , on the other hand, shows that the GT 
strength is almost independent of the tensor coupling, 
and is quenched by the same amount as in the case of 
the pseudovector coupling. 

When we take into account both the pseudovector and 
tensor coupling, the relativistic effect on the excitation 
energy of the GT state depends on their ratio, as in 
Eq. (|44|l , while the quenching of the GT strength remains 
in the same way as in other cases. This fact is shown by 
the factor (1 - 2v F /3) in Eq.®- 

In conclusion, the pseudovector and tensor coupling 
play a role to determine the excitation energy of the GT 
state in a similar way, but do not change the previous 
conclusion^, H that the GT strength in the relativistic 
model is quenched by 12% in nuclear matter and by 6% 
in finite nuclei, compared with the IFF sum rule value. 
The quenched amount may be increased by the cou- 
pling of the particle-hole states with nucleon-antinucleon 
states through the pseudovector coupling, as discussed at 
the end of the preceding section. Thepresent result is 
consistent with the recent experiment 12j, but in future 
work, other spin-dependent structure of nuclei, like the 
pion condensation and response functions of the charge- 
exchange reaction pi), should be understood consistently. 
It is also important to study the present relativistic ef- 
fects on the (5 decay for its precise discussions in neutrino 
physics. 
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